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Abstract 

The formation of dust ion-acoustic shocks (DIASs) in a four-component quantum plasma whose 
constituents are electrons, both positive and negative ions and immobile charged dust grains, is 
studied. The effects of both the dissipation due to kinematic viscosity and the dispersion caused 
by the charge separation as well as the quantum tunneling due to the Bohm potential are taken 
into account. The propagation of small but finite amplitude dust ion-acoustic waves (DIAWs) is 
governed by the Korteweg-de Vries-Burger (KdVB) equation which exhibits both oscillatory and 
monotonic shocks depending not only on the viscosity parameters ri± = fi±uj,p^/(?^ (where ^± are 
the coefficients of kinematic viscosity, ujp- is the plasma frequency for negative ions and Cg is 
the ion-sound speed) but also on the quantum parameter H (the ratio of the electron plasmon 
to the electron Fermi energy) and the positive to negative ion density ratio (5. Large amplitude 
stationary shocks are recovered for a Mach number (M) exceeding its critical value (Mc). Unlike the 
small amplitude shocks, quite a smaller value of r/+, r/_, H and [3 may lead to the large amplitude 
monotonic shock strucutres. The results could be of importance in astrophysical and laser produced 
plasmas. 
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1. Introduction 

Physically, a pair-ion plasma is similar to an electron-positron plasma, in which particles 
have the same mass and opposite charges. Negative ions are found to be an extra component 
(which may occur naturally or may be injected from external sources) in most space and lab- 
oratory plasmas [1,2]. There are a number of works for investigating such pair ion plasmas. 
Recently, Kim and Merlino [3] have discussed the conditions under which dust, when injected 
into a laboratory negative ion plasma, becomes positively charged for very large values of 
negative ion density > 500 times the electron density. Rapp et al [4] have discussed the 
possible role of negative ions in explaining their observations of positively charged nanopar- 
ticles in the mesosphcrc under nighttime conditions. Cooney et al [5] have investigated a 
two-dimensional soliton in a pair-ion plasma. Also, the role of negative ions in a laboratory 
dusty plasma have been discussed by Klumov et al [6]. An experimental investigation of 
the effects of negative ions on shock formation in a coUisional Q-machine plasma has been 
made by Luo et al [7]. Moreover, It has been pointed out that such pair-ion plasmas have 
potential applications in the atmosphere of D-region of the Earth's ionosphere, the Earth's 
mesosphere, the solar atmosphere as well as in the microelectronics plasma processing re- 
actors [8]. Takeuchi et al [9] in their work have reported the experimental observations of 
ion-acoustic shocks in an unmagnetized plasmas whose constituents are electrons as well as 
positive and negative ions. They observed wave steepening of positive or negative jumps 
for certain values of the positive to negative ion density ratio {f3 = n+o/n_o)- They sug- 
gested that their experimental observations can be well explained by means of the Korteweg 
de-Vries (KdV) equation in which the nonliner coefficient {A) can change its sign when a 
significant fraction of negative ions is present in an electron-positive ion plasma. Various 
laboratory experiments [10] have been conducted over the last few years to study the forma- 
tion of dust ion-acoustic shocks (DIASs) in dusty plasmas. Dust ion acoustic compressional 
pulses have been observed to steepen as they travel through a plasma containing negatively 
charged dust grains. Theoretical models [11] have been proposed to explain the formation 
of small amplitude DIASs in terms of the Korteweg-de Vries-Burger (KdVB) equation, in 
which the dissipative terms comes from the dust charge perturbations [12] and kinematic 
viscosity [13]. 

Quantum plasmas, where the finite width of the electron wave functions gives rise to 
collective effects [14,15], are currently an emerging field of research area. In view of its 
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potential applications in micro-electronic devices [16], nanoscale systems [17], in laser fusion 
plasmas [18 ], next generation high intensity light sources [19,20] as well as in dense astro- 
physical environments [21], various collective processes have been investigated in a number 
of research works [e.g. see Refs. 22-26]. 

There has also been much interests in investigating the strucure and dynamics of shocks 
in quantum like systems, such as nonlinear optical fibers and Bose- Einstein condensates [27- 
29] . The structure of such shocks are quite different from the classical ones where the shocks 
are typically governed by the transport processes, i.e., the viscosity and thermal conduction. 
Unlike classical fluids, quantum plasmas typically exhibit dispersion due to the quantum 
tunneling associated with the Bohm potential instead of dissipation [22-26]. For this reason, 
even a quantum shock propagating with constant velocity in a uniform media does not 
exhibit a stationary structure. Transition from initial to compressed quantum media occurs 
in the form of a train of solitons propagating with different velocities and with different 
amplitudes [22-26]. Such a train of solitons also provides a non-monotonic transition from 
the initial to flnal state of the medium. However, there are various quantum plasma systems 
in which both dissipation and dispersion play roles. Such important roles have been studied 
in different quantum plasma systems for the formation of ion-acoustic shocks (where the 
dissipation is due to kinematic viscosity) recently by Sahu et al [30] and Misra et al [26]. 

It is therefore of interest to examine the effects of kinematic viscosity as well as the 
quantum mechanical effects on the propagation of dust ion-acoustic waves (DIAWs) and to 
show how the dispersion caused by the charge separation as well as the density correlation 
due to quantum fluctuation, and the dissipation due to kinematic viscosity play crucial 
roles in the formation of shock waves instead of solitary wave solutions. We also investigate 
how the presence of negative ions modify the wave structures in a multi-component dusty 
quantum plasma. This is the aim of the present investigation. By using the standard 
reductive perturbation technique (RPT), the small amphtude DIAWs is described by 
the Korteweg-de Vries-Berger (KdVB) equation, where the Burger term appears due to 
kinematic viscosity determined by both positive and negative ions. The equation is then 
numerically solved to show that either oscillatory (dispersion-dominant case ) or monotonic 
(dissipation-dominant case) shock wave solutions are possible to exist depending on the 
nondimensional quantum diffraction parameter i7, the viscosity parameter 77-1- and the 
equilibrium positive to negative ion density ratio. 
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We have also recovered the large amplitude shock solutions for values of the Mach number 
(M) exceeding critical value (Mc) . 



2. Basic equations and small amplitude shock solutions 

We consider the propagation of DIAWs in an unmagnetized coUisionless quantum plasma 
composed of electrons, positive and negative ions and immobile negatively charged dust 
grains. The dynamics of DIAWs in our quantum dusty plasma is governed by the following 
set of hydrodynamic equations: 

e 90 I dpe d f I 92 



^^^^ — " v-e . 

me dx me^e dx 2me^ dx \ .Jnl dx^ ' ' 



dt dx 

~dt 



du+ Z+e d(p d'^u+ 



(3) 



du- du- Z^ed(j) d'^u- , , 

dt ^~ dx m_ dx ^~ dx'^ 

1^ = ^ (ne - Z+n+ + , (6) 

where riajUajfrLa arc respectively the density (with equilibrium value riao), velocity and 
mass for electrons {a = e), positive ions {a = +) and negative ions {a = — ); is the 
Planck's constant divided by 27r; 4> is the electrostatic wave potential; Pe is the electron 
pressure; x and t are respectively the space and time variables, and ii± is the coefficient 
of kinematic viscosity due to positive (negative) ions. At equilibrium the overall charge 
neutrality condition reads 

Ueo + Z_n^o + Zdo+rido = Z+n+o, (7) 

where Z±, Zdo are the charge states for positive (negative) ions and dusts, Udo is the equi- 
librium dust number density. We assume that the ions are cold, and electrons obey the 
following pressure law [31]. 
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where vpe — ^/S^bTVI/^ is the electron Fermi thermal speed, Tpe is the particle Fermi 
temperature given by ksTpe — /i^(37r^)^/^ngQ^/2me, ks is the Boltzmann's constant. Now- 
introducing the following normalizations 

X^CUp-X/Cs, t^Up-t, na^na/riao, Ua^Ua/Cs, (f) ^ e(f)/ (kBTpe) , (9) 

where a — e, +, — and Upa — y^naoe'^ / Soma is the a— particle plasma frequency, Cg — 



^/kBTpe/n^- is the quantum ion-acoustic speed, we get the following the normalized set of 
basic equations as : 

dx ^ dx 2/x dx \^/n^ dx'^ J ' 
dt dx ^ dx dx'^ 

(-) 

du- du- d(j) d'^u^ . . 

dt ^~ dx dx ^ dx'^ 

1^ = //rie - + (15) 

where = neo/ Z_n-o, j3 = Z_,_n+o/^-n_o connected through the charge neutrality condition 
[Eq.(7)] ji = l3 — \ — 5 with 5 = Zdondo/Z_n_o, m = m_/m^, rj^ = fi^-Up^/cj. and the 
nondimensional quantum parameter H = hvpe/ {kBTpf.) (the ratio between the electron 
plasmon energy and the electron Fermi energy) proportional to quantum diffraction. Now, 
integrating once the Eq. (10) with the boundary conditions viz. — > 1, ^ — > and — > 
at ±oo we have 

In order to investigate the propagation of small but finite amplitude DIAWs and to derive 
the required governing equation in our quantum dusty pair-ion plasma, we stretch the in- 
dependent variables as ^ = e^/^(x — \t),T — e'l'^t with r\± — €^^'^r]±o, while 77^0 is a finite 
quantity of the order of unity, while the dependent variables are expanded as 



na = 1 + enW + e^n^^^ + e^nf + • • • , 



(17) 



= + eui^) + e2«(f) + e^ii^f) + • • • , (18) 



= + £0(1) + eV^'^ + + • • • , (19) 

where o; = e, +, — and e is a small nonzero parameter proportional to the amplitude of the 
perturbation. Now, substituting the expressions from Eqs.(17)-(19) into the Eqs.(10)-(16) 
and collecting the terms in different powers of e, we obtain in the lowest order of e the 
dispersion law: 



^^^/2(l^ (20) 

Which shows that the DIAW can propagate outward or inward depending on the positive or 
negative sign of A, which increases (decreases) with (3 (/i). Now, in the next higher order of 
e we eliminate the second order perturbed quantities from a set of equations to obtain the 
required KdVB equation for DIAWs in our quantum plasma. 

dn , On ^d^n ^d'^n ,„ , 

where n = n^e^ and the coefficients of the nonlinear (^4), dispersive (B) and dissipative (C) 
terms are 



^ '1+^ _ 3^,(1 -m^P) ^ ^ 1^(16 - H% (22) 



77_o + mPrj+o 



2(l + m/3) 

Wc find that in Eq.(21) all the coefficients A , B and C arc modified by the inclusion of 
negative ions, while B and C are further modified by the effects of quantum diffraction and 
the effects of kinematic viscosity respectively. In order that the DIAWs propagate with finite 
velocity A we must have /x(= /3 — 1 — 5)>0or/3>l + 5. Note that in absence of viscosity 
term, Eq.(21) reduces to a usual KdV equation for the propagation of DIAWs, whereas for 
H — A (for which S = 0) it reduces to a purely Burger equation. 
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In order to investigate the nonlinear dynamics of DIAWs we have numerically solved the 
KdVB equation (21) for different sets of parameters. The equation (21) has a solution that 
can be represented as an oscillatory shock. However, when the dissipation overwhelms the 
DIAW dispersion and when the dissipative effect (kinematic viscosity) is in nice balance 
with the nonlinearity arising from the nonlinear mode coupling of finite amplitude DIAWs 
in a dusty pair-ion plasma, we indeed have the possibility of monotonic compressive or 
rarefactive shocks. The dynamics of the latter is governed by a Burger equation. In the 
numerical procedure the KdVB equation was advanced in time with a standard fourth-order 
Runge-Kutta scheme with a time step of 10~^ s. The spatial derivatives were approximated 
with centered second-order difference approximations with a spatial grid spacing of 0.2 m. 
The numerical solution of Eq.(21) is shown in Figs.l and 2 for different values of the system 
parameters /3, 77±, S and H. Figure 1 shows the oscillatory shock profiles of the electron 
density perturbation at the end of r = 0.3 s for m = 4, /3 = 4, 77+0 — O.3,?7_o = 0.2, 
6 — OA and H — 0.2. Even for small but non-zero values of H and r]±o we observe a train of 
oscillations ahead of the shock front decaying at ^ — > — 00. Increasing the role of quantum 
effects (H = 3) we observe a less number of oscillations ahead of the shock wave. The 
amplitude of the shock front seems to increase with the increasing values of H and f3, but 
decreases with S. The oscillations decay quite fast in case of relatively higher values of H and 
P and quite slow at higher S. Moreover, in case of large quantum effects {H — 3) and large 
dissipative effects (77+0 = 1-0, 77-0 = 0.9) and P — 10,S — 0.4 we recover monotonic transition 
from the oscillatory shocks (Fig. 2). For H < 3, we again find few oscillations ahead of the 
shock in which first few oscillations at the wave front will be close to solitons.The stationary 
solution of Eq.(21) can also be obtained analytically. We find that C always negative, B ^ 
according as H ^ 4. Since H > A corresponds to a lower density region (as can be found 
from the Fermi temperature-density relation mentioned earlier) we consider H < A, so that 
B > 0. Also, since m, (3 > l,A > 0. Thus, in our purpose A > 0,B > and C < 0. 
In order to find a stationary solution of Eq.(21) we use the transformation C — ^ ~ Uqt in 
Eq.(21), where Uq is the normahzed velocity of the DIA shock waves and obtain the following 
equation 

+ + - Uon + {Uo - ^A) = 0, (23) 
where we have imposed the boundary conditions n ^ 1, dn/dC,, d?"n/dC,'^ ^ as C — ^ 00. Eq. 
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(23) describes a shock wave [32] whose velocity in the moving frame of reference is Uq. Since 
the nature of the shocks depends on the system parameters (3, 'r)±, 5 and H, we consider the 
case where the dissipation term dominates over the dispersive term. In that case, Eq.(23) 
reduces to 

(.o-^^^cf. (24) 
Which yields upon integration the following monotonic compressive shock solution 

n=^[l-tanh{-^(e-t/or)}] (25) 

with the shock speed C/o,the shock height Uq/A and the shock thickness —C/Uq. Since, A 
increases with /?, and — C increases with /3 and 77^0, the shock height decreases as P increases, 
and the thickness increases with increasing values of l3 and ri±Q. When the dissipative effects 
are small, the shock will have an oscillatory profile, whereas for large values of r]±Q ~ 1, the 
shock will have a monotonic behavior as is seen in the numerical solution of Eq.(21). To 
determine the values of 77^0 a-nd (3 analytically (since C depends on both 77-1-0 a-nd (3) corre- 
sponding to monotonic or oscillatory shock profiles, we investigate the asymptotic behaviors 
of the solution of Eq.(23) for C ^ -oo.We substitute n{C) = 1 + iV(C), << 1 in Eq.(23) 
and linearize to obtain 

S^ + C^^ + (A-t/o)Ar = 0, (26) 
The solution of Eq.(26) are proportional to exp{p(), where 



It turns out that the DIA shock wave has a monotonic or oscillatory profile according as 
^ 2{A — t/o)-B.The stationary oscillatory solution of Eq.(21) is obtained as 



CC\ j: A-Uo 



7^ = 1 + Dexp I cos Ca/ — ^ | , (28) 



where D is a. constant and C = C ~ Uqt. 

3. Large amplitude shock solutions 
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We now consider large amplitude planar stationary shock. In the moving frame of refer- 
ence q — X — Mt, the basic normalized equations (10)-(15) may be integrated as 



1 d\/n: 



^^^^ 2/1 JrTe dq'^ 



(29) 



M 



M -u± 



(30) 



1 9 / H X , du+ 
Mu± + -ul = {-11, l)(t> + ?7±^ 



(31) 



= /xrie — + ri- 



(32) 



Eliminating the variables rie and and using the quasi-neutrality condition //rie = + 
n_ we obtain the following desired set of equations 



where 
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We have numerically solved the system of equations (33)- (35) by Runge-Kutta scheme start- 
ing from the initial conditions n+ — n- — 1.01, — 0.001 with a step size A^^ = 0.001. 
It is found that the perturbations develop into shock waves provided the Mach speed (M) 
exceeds its critical value = 5. For m = 4, j3 = 4,6 = 0.1, r^+o = 0.02, ?7_o = 0-01 and 
H = 0.4, the oscillatory shock solution is shown to exist for M = 6 in Fig. 3. As the value 
of H decreases, the number of oscillations ahead of the shock decreases together with the 
increase of the amplitude at the upstream side. The effects of 6 is to increase both the am- 
plitude and width of the oscillatory shock profile. Also, increasing the value of /3(say, 10), 
we can clearly see a less number of peaks and troughs at the upstream side = 0). Unlike 
the small amplitude case, the large amplitude DIA monotonic shocks can exist even for 
comparatively small values of ri±Q. As for example, for M = 6, m = 4, /5 = 4, 5 = 0.1, 
77+0 = O.2,77_o = 0.1 and H — 0.4 the monotonic shock profile is shown in Fig.4. Further 
increasing the value of the quantum parameter H, e.g., ior H — 1, the monotonic profile 
again transits into the oscillatory one (Fig. 5). We would like to stress that we have obtained 
both small and large amplitude stationary and nonstationary DIA shock structures in our 
quantum plasma model; the train of oscillations propagate along with the shock with the 
same velocity. 

4. Discussions and Conclusions 

In this paper, we have investigated the nonlinear propagation of DIAWs in a four- 
component quantum plasma composed of electrons, positive and negative ions and stationary 
negatively charged dust grains. Both the dissipative (due to kinematic viscosity) and dis- 
persive (due to Bohm potential) effects are taken into consideration for the formation of 
DIA shock structures. Both the stationary and nonstationary shock solutions are recovered 
in our quantum plasma model. The small amplitude nonstationary and stationary DIA 
shock structures are obtained numerically and analytically from the KdVB equation. The 
transition from oscillatory to monotonic shocks strongly depends not only on the quantum 
parameter if and the viscosity parameter r]±, but also on the positive to negative ion density 
ratio P (which is mainly due to the inclusion of negative ions) . The shock structures are also 
modified by the effects of H, j3 and 5. Numerical simulation reveals that the small amplitude 
nonstationary shock (monotonic) transition occurs compratively at larger values of H{> 3), 
/?(> 8) and r]± ~ 1, whereas for the large amplitude shocks comparatively lower values of 
if(~ 0.4), /3(~ 4) and 77±(~ 0.2) are required. Such significant modifications of the shock 
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wave front structures in our quantum plasma model could be of interest in astrophysical 
and laser produced plasmas. Furthermore, it is suggested that the experiments should be 
designed to look for dust ion-acoustic shocks (DIASs) and solitons in such quantum dusty 
pair-ion plasmas. 
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Figure captions 

Figure 1: Small amplitude oscillatory shock profile [nonstationary solution of Eq.(21)] for 
m = 4, /3 = 4, r/+o = 0.3, r/_o = 0.2, S = 0.4 and H = 0.2. 

Figure 2: Monotonic shock structure [nonstationary solution of Eq.(21)] for m = 4, 
(3 = 10, r)+o = 1-0, r)-o = 0.9, S = 0.4 and H = 3. 

Figure 3: Large amplitude oscillatory shock profiles for M = 6, m = 4, /3 = 4, 5 = 0.1, 
77+0 = 0.02, 77_o = 0.01 and H = 0.4. 

Figure 4: Large amplitude monotonic shock transition from the oscillatory profile in Fig. 3 
for M = 6, m = 4, /3 = 4, 5 = 0.1, 77+0 = 0.2, r]_Q = 0.1 and H = 0.4. 

Figure 5: The monotonic shock profile in Fig. 4 transits into the oscillatory one for 
H = 1. Other parameter values remain the same as Fig. 4. 
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